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Abstract 

Gyrokinetic field theory is addressed in the context of a general Hamiltonian. The background 
magnetic geometry is static and axisymmetric, and all dependence of the Lagrangian upon dy- 
namical variables is in the Hamiltonian or in free field terms. Equations for the fields are given 
by functional derivatives. The symmetry through the Hamiltonian with time and toroidal angle 
invariance of the geometry lead to energy and toroidal momentum conservation. In various levels 
of ordering against fluctuation amplitude, energetic consistency is exact. The role of this in under- 
pinning of conservation laws is emphasised. Local transport equations for the vorticity, toroidal 
momentum, and energy are derived. In particular, the momentum equation is shown for any form 
of Hamiltonian to be well behaved and to relax to its magnetohydrodynamic (MHD) form when 
long wavelength approximations are taken in the Hamiltonian. Several currently used forms, those 
which form the basis of most global simulations, are shown to be well defined within the gyrokinetic 
field theory and energetic consistency. 
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I. INTRODUCTION 



Gyrokinetic theory is a well founded formalism by which particle motion is treated in 
terms of drifts of particle gyrocenters rather than the combination of gyromotion and drifts 
of the particles. Particle motion in a magnetic field is set up with a drift kinetic Lagrangian 
assuming arbitrarily large magnetic field scale length [l-3|, and then Lie transform tech- 
niques assuming small product of gyroradius and field amplitude are applied to obtain a 
Lagrangian still independent of gyrophase angle but valid for a gyroradius, while still small 
compared to background scale lengths, of arbitrary order with respect to the scale of E- 
cross-B eddies j4-7|. Alternatively, the field variable amplitude may be left arbitrary when 
the small parameter for expansion is the local ratio between the gyroradius and the scale 
of potential variation as. Initially in this formulation, a back transformation was used 
to obtain the self consistent field polarisation equation for the electrostatic potential js], Q]. 

brm are an application of 



i 
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More recently 



Taken together the gyrokinetic Lie transform and this back trans 
push forward and pull back transforms in differential geometry 
the entire Lagrangian is set up as the integral of a Lagrangian density over the phase space, 
with the polarisation equation obtained as the Euler-Lagrange equation by varying the elec- 



trostatic potential in the field Lagrangian 
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13| . These two approaches were shown to 



11| . Moreover, in an analysis 



be equivalent in the recent review by Brizard and Hahm 
of gyrokinetic transformation of the general Landau collision operator, the method of Lie 
transforming the Lagrangian and deriving the Euler-Lagrange equations was shown to be 
equivalent to the Poisson bracket transform of the Vlasov or Boltzmann kinetic equation di- 
rectly, with the latter method able to treat coUisional dissipation but with the Lie transform 
useful in deciding which coordinate map to use 1^. Following demonstrative gyrokinetic 
simulation of the internal kink mode 15|, gyrokinetic theory was explicitly linked back to 



MHD for long wavelength electromagnetic oscillations and instabilities |16l-ll9l|. Using the 
large amplitude/long wavelength form of Ref. (9| the correspondence to nonlinear reduced 
magnetohydrodynamics (MHD) and with the long wavelength limit of the small amplitude 
version was shown at the level of the Lagrangians. The theory is now a fully self consistent 
Lagrangian field theory. As it does not depend on assumptions concerning the form of the 
distribution function, it is also necessarily a total-f formulation. 

Energy conservation has been well known since the development of the Lie transform 
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version of the theory as cited above. Momentum has received less attention despite the 



general demonstration of conservation via the Noether theorem 
rise of gyrokinetic treatments of neoclassical flows 
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12|. However, since the 



discussion of the Coriolis drift 



and turbulent equipartition effects 22|-|25| has emerged during observation of the tokamak 
momentum pinch 26|. The turbulent equipartition scenario as an indirect effect involving 
exchange channels more than actual drive effects especially underscores the role of strict 
conservation in a complicated, nonlinear physical situation. Momentum conservation was 
also demonstrated for evolution of axisymmetric flows and currents toward equilibrium in 



the context of total-f electromagnetic theory and computation under edge conditions [27] . 

Gyrokinetic field theory was not necessary to build the original gyrokinetic computational 
models whose self consistent equation for the field potential was referred to as the gyrokinetic 
Poisson equation 28l-l30l|. Despite recovery of these equations by the Lie transform and 



field theory methods cited above, the usefulness of the gyrokinetic Poisson equation for 
determining the electric field has been criticised 3l|. Within the field theory, however, the 
Euler-Lagrange equation for the electrostatic potential is indeed one and the same with this 
gyrokinetic Poisson equation 12|. Obtaining the equation via different methods outside of 
the field theory would appear to break the inherent consistency unless the method were 
found to be equivalent to use of the field theory, as the earlier versions in fact are. 

We therefore address the general question of momentum conservation within the energetic 
consistency afforded by the field theory version of gyrokinetics. The results of any particular 
Lie transform are assumed to be given, with the only stipulation being the general form of 
a Lagrangian in which all dependence on time dependent field variables is transformed into 
the Hamiltonian {i.e., the "interaction Lagrangian" as described in basic texts [32j). Both 
energy and momentum are considered, with the role of symmetries as paramount. In a 
tokamak under low frequency conditions, the conserved momentum is toroidal momentum; 
the other components of the momentum vector use the spatially dependent background 
magnetic field as an anchor. We leave the issue of a conserved momentum vector, a tensor 
transport flux, and the anchor to future work and concentrate on toroidal momentum. The 
result is that momentum transport and conservation yields familiar content and transport 
fluxes, generally and even in the specific cases of conventional models. The route back to 
MHD is explicitly shown via appropriate choice of the Lagrangian/Hamiltonian. The main 
results are valid for any ordering scheme which might be used since their demonstration 
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does not depend on ordering but uses general functional form of the Hamiltonian on the 
field variables and a field term describing shear- Alfven disturbance magnetic energy. 

Following sections of the paper describe the field theory model in general with emphasis 
on energetic consistency, recover global conservation laws via the Noether theorem, address 
the conservation of generalised vorticity, and then address toroidal momentum and energy. 
A 4-dimensional antisymmetric bracket form of the gyrokinetic equation is derived (the 5th 
and 6th dimensions do not enter due to the conservation of generalised magnetic moment 
and the lack of dependence by any dependent variable on gyrophase angle), which not only 
greatly facilitates the mathematics but is also suitable for computations 27(]. Two comments 
on orderings are given, both how naive ordering schemes can violate energetic consistency, 
and on a result showing that momentum evolution via small fluctuations occur in two places 
via terms of the same order {i.e., orderings that occur there enter order by order in the same 
term multiplied by factors of order unity). The general momentum and energy transport 
equations that are derived are independent of any ordering. The route back to MHD is 
shown using what can be called an "MHD Hamiltonian" is shown explicitly for momentum 
and a mean field fluctuation model within that is given. The salient mathematics regards 
application of functional derivatives, for which the background and main operations are 
given in Appendix A. Then, Appendix B treats conventional models and their equivalence 
to gyrokinetic field theory versions in each case via appropriate choice of the Lagrangian. 
The results on momentum depend on a cancellation mandated by the transport equation 
for generalised vorticity; Appendix C gives the version of this by considering the torque due 
to a charge source under quasineutral conditions. 

II. DESCRIPTION OF THE GYROKINETIC FIELD THEORY MODEL 

The general gyrokinetic theory follows from Lie transforms applied to the extended phase 
space Lagrangian for the particle gyrocenter motion {4-^. A phase space kernel locating 



the particles transforms the Lagrangian to a Lagrangian density [12|]. In the electrostatic 
version the phase space integral over this density is the entire action. In the case of an 
electromagnetic model the magnetic energy is added as a free field contribution as in general 
electrodynamics (for background see the text by Landau and Lifshitz [32]). The electric field 
energy contribution is neglected since in a magnetised plasma the ExB kinetic energy of the 
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particle drifts is much larger; this is the same statement as quasineutrality since a zero 



space charge density is the natural result 



12| . Variation of this action with respect to 



the gyrocenter coordinates gives the Euler-Lagrange equations for the particles. Liouville's 
theorem is used to convert these into an equation for a distribution function, which serves as 
the gyrokinetic Vlasov equation. Variation of the action with respect to the field potentials 
gives the polarisation and induction equations giving the self consistent response of the 
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13|. 



electrostatic potential and parallel magnetic potential, respectively 

Gyrokinetic theory was recently given in terms of the field theory for transonic ExB flows 

9|, readdressing the large amplitude (strong ExB flow) version of Ref. jsj. Following Ref. 

33| . the ExB flow resulting from the appearance of the large-scale potential in the lowest- 
order Euler-Lagrange equations was used in the coordinate transformations, rather than a 
background flow given by a fluid analysis. In our case, however, the same potential was used 
for all flow dynamics without splitting the fluctuations from the background; the dependent 
field variable was used for both purposes. The choice of Lie transform was then changed 
to move all effects from the flow potential into the Hamiltonian. The model was shown to 



35[, for weaker 



recover the conventional ones [6|, LZj, including nonlinear reduced MHD [34 
flow amplitude and larger scale. 

The particle Lagrangian itself starts with the form from electrodynamics and is Lie trans- 
formed into a low-frequency form in an expansion which formally uses the amplitude of the 
drifts as a small parameter; this can be either the fluctuation amplitude as in Refs. js, 6| or 
the gyroradius compared to the dynamical scale length as in Refs. ^, 9|. 

The Lie transform method is rather general, and various choices exist, but for present 
purposes it is useful to know that the choice can always be made to arrange the particle 
Lagrangian, Lp, so that the symplectic part (the vector of coefficients, pi, of particle coor- 
dinate time derivatives, q\ in the representation Lp = picf" — H) depends on background 
geometry only, and all time dependent field effects appear in the Hamiltonian, H, only. 
Since the latter serves as the time component of the underlying fundamental one-form, this 
casts the Euler-Lagrange equations in a form where partial time derivatives on dynamical 
field variables are absent. Resulting geometric quantities, including volume elements and 
Jacobians, are strictly static. This was the Lie-transform strategy used in Ref. [9[], empha- 
sising the general desirability of arranging Lp with all time dependence in H not only for 
computations but also in proving correspondence to other forms such as nonlinear reduced 
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MHD. For the case of tokamak geometry, this has the added benefit of isolating all toroidal 
angle dependence into H as well, which will be seen to facilitate the proof of results involving 
toroidal momentum conservation. 

We therefore assume a particle Lagrangian which for any degree of expansion has been 
Lie transformed into the following form 



where {Zp} = {R,^^, /i} are the particle coordinates (gyrocenter position, parallel canonical 
momentum, gyration magnetic moment), ?9 is the ignorablc sixth coordinate giving the 
gyrophase angle, e and m are the species charge and mass, A and b are the potential and 
unit vector for B = VxA the background magnetic field, and H is the Hamiltonian. The 
canonical parallel momentum version is used, so that the time dependent parallel magnetic 
potential A\\ appears only in H. All flow dynamics due to the time dependent electrostatic 
potential (f) appear also only in H. Here and below, the parallel subscript denotes the 
component locally parallel to B. 

The Hamiltonian depends on both field potentials and A\\, evaluated at the gyrocenter 
positions via, e.g., 0(R), as well as the gyrocenter phase space coordinates. 



The parallel velocity U is not used explicitly as a coordinate but can be defined as a derivative 



It is important in all the derivations to note that U has spatial and time derivatives through 
its dependence on 

The dependence of H upon the fields 0, Ay involves differential operators such as spatial 
derivatives which commute generally with variations or spatial or time derivatives, or the 
gyroaveraging operator Jq which depends on /i and B as well as spatial derivatives. One 
has to know whether these operators commute with derivatives. With all time and toroidal 
angle dependence transformed into H we ensure to be working with a representation in 
which this is true generally for differentiation with respect to either time or toroidal angle 
although not for all components of the gradient operator. 




(1) 



H = H{R,p„ii,(P,A\\) 



(2) 



oiH, 




(3) 
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Formally, Jo has the form in wavenumber space of multiplication of Fourier coeffi- 
cients by the zeroth Bessel function Jo{k±pL), where pi is the particle gyroradius given 
by Pl = v±/\eB/mc\, with gyrofrequency \eB/mc\, or in terms of the coordinates by 
p| = 2pB /[m{eB/mc)'^]. The perp subscript denotes the component in the plane locally 
perpendicular to B. Hence Jq is time symmetric in any geometry but toroidal angle sym- 
metric only in tokamak geometry. The Jq operator may be cast as a series of perpendicular 
Laplacians, V^^, so in the local transport equations to be derived below it is sufficient to con- 
sider H with arbitrary dependence on field amplitude and the field gradient and Laplacian. 
This becomes necessary when considering the role of functional derivatives in the theory. 

The particle equations of motion are found from the Euler-Lagrange equations resulting 
from Lp. The drift motion, 

B*.— = VH--- + UB* B;^ = -B*-VH (4) 

separates naturally from the gyromotion, 

dp, ^ dd e OH , , 

dt dt mc dp 

where some standard definitions are 

A* = A + p^-b B* = VxA* Bl = hW (6) 

Drift tensor notation is used, with 

F = VA - (VA)'^ (7) 
where superscript T denotes the transpose. It follows that 

F = e-B Vxb = -V-^ B* = B-j9,V--^ (8) 

B e B 

where e is the rank-three Levi-Civita pseudotensor. Phase space volume conservation is 
expressed by 

5|;-(v^S;Zp)=0 (9) 

where ^/gis the determinant of covariant components of the coordinate metric. With V-B* = 
0, Eq. (E]) implies 

^ cF dB* 
7 



Guaranteed by the definition of B*, tliis also determines tliat tlie quantity -By serves as tlie 
volume element of the velocity space, where the volume element of the entire phase space is 

This particle Lagrangian is converted to a system one by placing the particles in phase 
space via the kernel G{Zp, Z) with 

Z -> {x, z, w} Zp {R, p,, /i} (11) 

giving the correspondence between phase space coordinates Z and gyrocenter coordinates 
Zp, respectively. Due to some notational difficulties with the rest of this work, however, we 
will dispense with this distinction between Z and Zp, leaving the role of G to be understood 
once we already have the particle equations of motion in Eqs. (HfSjl . 

The phase space integral is denoted / dA. The integration domain dA is given as a com- 
bination of the velocity space and configuration space domains. These are given respectively 
by 

dA = dV0 dW dV = ^dx^ dx^ dx^ dW = 27rm-^ dp, dfi 5[f (12) 

where ^/g is the determinant of covariant components of the coordinate metric, and not- 
ing that the form of dW is determined by phase space conservation (the i?y factor) and 
normalisation (the 27rm~^ factor). The Lagrangian for the entire particles/field system is 
then 

i = E/<'A/L.-/^V§ (13) 

bp 

where the sum is over species. The electrodynamic field term {E'^ — B'^)/87t reduces to 
—Bj_/87r as the assumption of quasineutrality eliminates E'^/Stc in favour of the ExB kinetic 
energy of the particle drifts, and the assumptions of low frequency u ^ k±VA and low plasma 
beta P = Snp/B"^ ^ 1 restrict magnetic variation to the parallel magnetic potential A\\. A 
simplified version usable for the present purposes is 

(14) 

defined with the perpendicular spatial derivatives. Hence Bj_/8TT is identified as the energy 
in shear Alfven magnetic disturbances perpendicular to the equilibrium magnetic field. More 
general forms are possible (cf. Appendix A) but the field term is always quadratic in A\\. 
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Under conditions of magnetic compressibility, enters as well (cf. Sec. Ill C of Ref. 
but the structure of the theory remains as presented herein. 

The gyrokinetic Vlasov equation is found via variation of the Zp components according 



to characteristic methods 



12| . or equivalently by application of Liouville's theorem to the 



particle motion in Eqs. yielding 

Derivatives with respect to /i or do not appear, because / and H are independent of 
and ^ is conserved in the gyromotion. Once we have the gyrokinetic equation in this 
form, the distinction between phase space and gyrocenter coordinates may be left implicitly 
understood, since we no longer consider particles or gyrocenters as discrete entities. 

The term "drifts" refers to the drift motion described by Eq. (jlj), especially the spatial 
part R. The part resulting from the field dependent variables is entirely contained in H. 
Hence when we refer to the treatment of drifts we mean the construction of H and in 
particular drifts to a certain order means the contributions to H due to an expansion up to 
that order. The results we will obtain do not depend on the form of ordering (just on the 
functional form of the dependence of H upon </> and ^n), but at certain points we will need 
to refer to the result in terms of a certain ordering. 

The equations for the fields are determined by functional derivatives (cf. Appendix A and 
the background references cited there). The self consistent polarisation equation (also called 
gyrokinetic Poisson equation) is given by the Euler-Lagrange equation for cf) from this same 
Lagrangian. It is found by varying the Lagrangian with respect to 0, yielding an integral 
over dV of 50(x) times a coefficient, which is required to vanish. It is the same statement 
as requiring the functional derivative of L with respect to (j) to vanish. This produces 

E^ = (16) 

The species-summed functional derivative of fH vanishes alone because appears only in H. 
The functional derivative implies velocity space integration because it is defined with respect 
to the space covered by dV. The functional derivative combination yields the gyrokinetic 
charge density for the particular H used (the assumption of quasineutrality sets this to zero). 

The self consistent induction equation (also called gyrokinetic Ampere equation) is ob- 
tained by variation of the field potential An. It is the same statement as requiring the 



functional derivative of L with respect to A\\ to vanish. This produces 

°p II 

with the field term appearing on the right side arising from the field term —B'j_/8TT in L. 
The functional derivative combination yields the gyrokinetic current (times — 1/c) for the 
particular H used. 

A. Antisymmetric Bracket Form of the Gyrokinetic Equation 

It has been found previously that maximal symmetry in the representation of th e gy roki- 
netic equation (Eq. [T5|) is helpful to the understanding of the conservation laws 27|. We 
observe that 

dA* c d cF 

^ = -h hence ^e-A*='-^ (18) 
opz e opz e B 

If we define 



G = e ■ A* (19) 



we may recast Eq. f|T5l) as 



This has the structure of one 3-bracket of indices {ahz} 

[H, fU = -g^[H, fU + C^aG^'m + {VtC^'m, fU (21) 

for each pair of spatial coordinates {ab} with index z denoting the pz coordinate. This can 
also be written as 

[H, fU = ^"'^ ( ^[^' Z]'^^ + (VaA:)[^, + (V,A:)[/J, f]za\ (22) 

where on the right side the Einstein summation convention is used for repeated (up/down) 
indices. In each case the 2-bracket form is 

[H,fU = H,aU-H,,f,a (23) 



with the comma denoting differentiation with respect to the coordinate whose index is given 
by the subscript. Since there is no p^-component of A* we may add 3 more fictitious 
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3-brackets, one for each pair of spatial indices and pz with A*. The entire combination 
becomes 

^ + S'^'^'HJ^.Al, = (24) 

where S is the rank-four Levi-Civita pseudotensor in the 4-space covered by dV ® dp^. The 
components of e"'"^ are times ±1 or depending on the permutation of spatial indices 

{abc}. The components of f!"^^*^ are l/^/gB*^ times ±1 or depending on the permutation of 
indices {abed} in the 4-space domain. The 3-space order is {123} for dx^ dx^ dx^ and hence 
the 4-space order is {1232;} for dx^ dx'^ dx^ dp^. Positive, negative, and zero permutations of 
these give the other components. 

It was previously observed that axisymmetric momentum conservation follows directly 
from this form of the gyrokinetic equation, simply due to symmetries in the indices [23]. In 
this work this antisymmetric bracket form will be used to facilitate proof of the conservation 
laws for energy and toroidal momentum for general dependence of the Hamiltonian upon 
the dynamical fields. 



III. GLOBAL CONSERVATION LAWS 



The conserved energy is found from the total action 
theorem, applying small variations to the time component 



(it L(Zp, 0, All, t) via Noether's 



32| . In the gyrokinetic case this 



12|, and 



has been done before, both from the discrete particle characteristics point of view 
from the continuum/field representation using constrained variations [l3|. Since L is first 
order in all the time derivatives this becomes the combination of all the picf" terms less the 
Lagrangian. This defines the Noether energy as 



S = J2 f dAfH+ [ dV 

sp 



Stt 



(25) 



In the electrostatic case it is simply the integral over fH summed over species. Since the 
background magnetic field (through A) is not varied, it does not appear in the Noether 
energy. 

The same follows for the Noether momentum, by applying small variations to the space 



component 32| . In the gyrokinetic case this was given in abstract form by Refs. 12|, 



13|, but 



not for specific cases. Working out the space components is complicated by the fact that 
in low frequency dynamics in a magnetised plasma, the magnetic field serves as an anchor 
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for momentum, so that the general four-vector version is not conserved for the dynamics 
under consideration: not only is the Poynting momentum neglected against the plasma 
momentum in ExB motion (through the assumption of quasineutrality) , but also the neglect 
of compressional Alfven dynamics removes the exchange with the background field. In a 
tokamak, only toroidal momentum is conserved. For L of the form given in Eq. (1131) . the 
conserved toroidal momentum is simply given by the toroidal canonical momentum weighted 
by / and summed over species. The Noether toroidal momentum is 



V = Y.jdKfP^ (26) 



sp 

where = dL/dip and (p is the geometric toroidal angle. This result is a consequence of 
all the dependence of L upon ip is in the time component [H] or in the field terms and in 
the latter there is no time derivative dependence. 

One other consideration is that one would like a local form of the conservation law in 
terms of a vector momentum density, a symmetric stress tensor for momentum transport, 
and a vector describing the magnetic field anchor explicitly, but this has yet to be worked 
out. Herein, we consider energy and toroidal momentum only, and explain their conserva- 
tion using the antisymmetric bracket form of the gyrokinetic equation and the functional 
derivatives which describe the self consistent field equations. One motivation for this is 
that it is possible to directly code the antisymmetric bracket form in numerical simulations, 
so it is then known that the form of the equations as actually used is indeed energetically 
consistent. 

Once the Noether energy and toroidal momentum are known, appropriate operations 
on the equations of motion (here, the gyrokinetic equation and the self consistent field 



equations) may be used to construct a local form with the time derivative o 



■ an evolving 



energy/momentum density and the divergence of an overall transport flux 13|. We will do 
this herein as well, as part of the overall motivation to establish the correspondence to fluid 
and MHD forms. 



A. Energetic Consistency 

Both of these approaches lead naturally to the known results on energetic consistency, 
namely that the same H must be used to obtain the gyrokinetic equation and the (related) 
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field equations. Approximations are done in L (hence H) and then the equations are derived 
without further approximation. Specifically, orderings in the derivation of L are used (via 
Lie transforms or some other method), but not thereafter in the derivation of the Euler- 



Lagrange equations. Ref. 



12j points out specifically that if the drift motion is to be followed 



with lowest-order forms of H, all higher-order forms must be cast into field terms. For 
example, ii H = Hq + Hi + H2, with the last piece containing the ExB energy, and it is 
desired to advance the gyrokinetic equation only with Hq + Hi, then the term fH2 in L must 
be replaced by /o-f^2 where /o is a background static form which can be thought of as part 
of the geometry. Then, since H2 does not multiply /, it is not involved in the gyrokinetic 
equation itself but only as a field term which would appear on the right hand side of the 
polarisation equation, e.g., 

T.i^fiHo + Hi) = -Y.yoH2 (27) 

This is referred to as linearised polarisation. The two assumptions go together: first order 
drift motion, and the appearance of /o in the polarisation term (the right side of Eq. 1271) 
Conversely, if one desires to keep the dependent variable / in this term, restoring Eq. (fT6|) . 
then the corresponding H2 must be kept in the drift motion. This is the basic statement of 
energetic consistency in a total-f global model and the essential references 
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13| arrived 



at this result ten years ago. The same result is found for the same reasons in gyrofiuid field 
theory models which have a different starting point but are also Lagrangian/Hamiltonian 
models {sg], [s^]. It is related to the connection in fluid models between advection and 
divergence forms of the equation of motion with respect to the polarisation drift velocity 
in a fluid model (why the polarisation drift must be kept in advection if the species mass 
density involves the dependent variable for species particle density 38|). 

A clear extension of this is that in any discussion of drift motion past first order, say to 
order n, the / must be kept as the dependent variable in all of the terms Hq + Hi + ■ ■ ■ + 
Hn in the functional derivatives in the polarisation equation. Any energetic contributions 
Hn+i + ■ ■ ■ must then either be dropped or combined with a background /o in L, with the 
polarisation equation then becoming the appropriately generalised version of Eq. fl27|) . Of 
course, if L is derived or constructed first and then the Euler-Lagrange equations are derived 
without approximation thereafter, then energetic consistency becomes a guaranteed result. 
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B. Time Symmetry and Energy Conservation 



Using the antisymmetric bracket form of the gyrokinetic equation (Eq. we multiply 
by H and use the linearity of the derivatives H f i, = {fH) h — fHh and the antisymmetry 
(the form with H aH f, vanishes due to the permutation of indices in S"''"^'^), to find 

^ifH) + S'^'^'HMH),Al, = (28) 

which is the local energy equation in phase space. 

Integration over phase space and summation over species yields 

El^4(m-i:j'"^ff (29) 

with the bracket vanishing under the integral. Under the integral the right side is replaced 
by functional derivatives (cf. Appendix A) 

E/->)^/-Ef^f./-E|^^ (30) 

The first term on the right side vanishes, due to polarisation (Eq. [T6|) . The second is replaced 
by the field term in A\\, due to induction (Eq. [T71) . so that 



sp 



Integration of the divergence operator in V]_ by parts then yields 
Identification with in Eq. (|T^ then recovers 

which is the same as obtained using the Noether theorem (cf. Eq. 1251) . This is the statement 
of energy conservation and it is valid under time symmetry for any dependence of H upon 
(f) and A\\ given the form of L stated above. 

The dependence of energy conservation upon time symmetry is contained in the step 
from the time derivative to the functional derivative, as d/dt must commute with any of the 
differential operators involved in the dependence of H upon and A\\. The requirement of 
energetic consistency is evident in the fact that the same H is used in the gyrokinetic equation 
as that whose functional derivatives appear in the polarisation and induction equations. 
Also, the same / must appear with each of the terms in H in all cases, or else the symmetry 
is broken. 



14 



C. Axisymmetry and Toroidal Momentum Conservation 



In this case the properties are different. For energy, H depends on both time and toroidal 
angle, but the symmetry of the bracket allowed combination of {fH) there. Canonical 
momentum at the particle level is given by 

P^^-^A*^^-^A^+p,h^ (34) 

This form is both static and axisymmetric, but does not appear in the bracket. We multiply 
by to find 

+ S^'^'HMP^hAla = f-/^'''H^A;^,Al, (35) 

Because the two appearances of A* in the right side appear with different indices, we may 

make some symmetry arguments. First, the indices b and d cannot be (/?, due to axisymmetry. 
Second, index c cannot be v?, or else A* appears twice and the remaining permutation over 
indices {abd} vanishes. Hence, index a must be (p and the others are among the other three 
coordinates, so that 

a f)]-f p 

Of the remaining terms, index b cannot be z because there is no A*, so in each term one of 
a or c is 2; while the others are the two coordinates covering the perpendicular plane, which 
we can label 1 and 2 (hence dx^ in dV is dip). We also observe that 

which cancels the (e/c) factor. Noting that the units of are 1/ .^/gBp we set the positive 
permutation as {(p zl 2} and permute the {z 1 2} indices to find 

ls^^'^A;^aAl,, = [b^{Al, - Al,) + b,{A*2,^ - a;,) + b,{A;^, - A^^)] (38) 

where we have eliminated the two (zero) terms A* ^^* 2 and A* 2^*,i' replaced them 
with the two (zero) terms &2^i,(^ and biA2^^, respectively. We observe that 

[b^{Al, - Al,) + bM;,2 - A*, J + hiAl^ - A;^,)] = b • Vx A* = b; (39) 

and note the switch in the order of coefficients to the expression just above. Therefore we 
have 

l£^''"'A;^Al,c = -1 (40) 
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Putting this into the right side of Eq. (136|) , we find 

f) f)J-f 

^(/P,) + S^'^%MP.),i>Al, = -f^ (41) 

which is the local toroidal momentum equation in phase space. 

This is the same form as the result obtained from a canonical representation of the particle 
Lagrangian 

TTIC 

Lp = P^^ + Pee + P^if + —fi^ - H (42) 

written directly in terms of the coordinates {ip 9 ip} as was once usual 3|). The corresponding 
Euler-Lagrange equation for the toroidal angle is 

dH 

dip 

due to the axisymmetry of the rest of Lp. Using the advection forms 



= -— (43) 



we find 

3^(W) + Zp-^(W = -/^ (45) 

which is the same form as in Eq. fl41l) with the bracket recast in terms of an advection. 
Returning to Eq. fHT]) . integration over phase space and summation over species yields 



S/<;A|(/P.)^-E/^A/f (46) 

with the bracket vanishing under the integral. Under the integral the right side is replaced 
by functional derivatives (cf. Appendix A) 



i:/rfA|(/P,)=/^VE^|^+/^VE§^§i (47) 



d_ 

di 



The manipulations follow the energy derivation, with d/dt replaced by 8/ dp. In this case 
the derivative of with respect to p vanishes under the integral, so that 

Y,fdAfP^ = (48) 

sp 

which is the same as obtained using the Noether theorem (cf. Eq. 1261) . This is the statement 
of toroidal momentum conservation and it is valid under axisymmetry for any dependence 
of H upon and A\\ given the form of L stated above. 
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The dependence of toroidal momentum conservation upon axisymmetry is contained in 
the step from the toroidal angle derivative to the functional derivative, as d/dip must com- 
mute with any of the differential operators involved in the dependence of H upon and A\\. 
The requirement of energetic consistency is evident in the same way as for energy with the 
same loss of consistency if the symmetry between functional derivatives and the gyrokinetic 
equation is broken. 



D. Phase Space Continuity Forms of Energy and Toroidal Momentum Conserva- 
tion 



Eqs. fl28|4ip give the antisymmetric bracket forms of the local energy and toroidal mo- 
mentum equations in phase space. Using toroidal momentum as an example, We identify 

e^'^'^H^j^.Ai, = K-^^ (49) 

Then, the phase space conservation condition (Eq. [9]) can be used to express Eq. (HTj) as 

which is the phase space continuity equation for toroidal momentum. 
Similarly, we may express Eq. ( I28l) as 

which is the phase space continuity equation for energy. 

Eq. ( 150|) makes it obvious that /P^ is conserved locally in axisymmetric systems {8/ dip = 
0), and globally in any geometry. In Eq. ( 15T|) the term appears instead with dH/dt, which 
eventually accounts for the magnetic energy as in Eq. ( l33l) . We will use these continuity 
equations to produce local transport equations for both momentum and energy later. 



IV. A COMMENT ON ORDERING 



Before the emergence of the Lagrangian/Hamiltonian approach to drift kinetic particle 
motion [l-^, in which the equation for / is built explicitly using the phase space positions 
of the gyrocenters as dependent variables, it was customary to start with the Vlasov (or 
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Boltzmann) equation for the particles themselves and apply gyroaveraging through a suc- 



cessive ordering 



39 



40|. This causes problems, however, if applied naively (and strictly) 



to the field equations. Considering an electrostatic model {A\\ = 0), we may expand H in 



terms of small amplitude fluctuations (e0/Te ~ = pJL 



gyroradius and L±_ is the profile gradient scale length 4l|, |42| ) 



where pi is the thermal ion 



H = H^ + Hi + H2 + --- + Hn (52) 

where at each order n the term if„ is n-th order in 0. As long as / is not expanded order by 
order, there is no problem. The same / multiplies each if„ in turn and due to the linearity 
property the functional derivatives add, producing Eq. f lT6|) term by term. 

However, if / is also ordered such that /o is the background (usually a Maxwellian) and 
/i is the fluctuation, then there is a problem. Recall that if drift motion is included to 
order n in H through then the polarisation equation must include in the functional 
derivatives to preserve energetic consistency. The / must include both /o and /i. However, 
if the ordering is truncated at order n then the term 6{fiHn)/S(j) is missing. Formally, it is 
order n+1. So this (n+l)-th order term must be kept, but in doing so we violate the ordering 
scheme. If the ordering scheme is applied to expand / and keep all terms up to order n, 
dropping all order n + 1 terms, then this one piece will be missing. This problem is present 
at any order of expansion, at the last order. The only acceptable solution for orderings 
is to expand H in orders but not /. That is, polarisation is not to be linearised (cf. Eq. 
|27|) if contributions above linear order in H to the drift motion are considered. Hence any 
discussion of orderings in which higher order drifts (even second order) are considered should 
be done under full energetic consistency. As noted, the field theory version of gyrokinetics 
is the only straightforward way to guarantee this. 



V. THE EXB VORTICITY TRANSPORT EQUATION 

We may form an equation for the gyrocenter charge density by multiplying the gyrokinetic 
equation by the charge e for each species and summing over species. For general dependence 
of H upon (f), the terms linear in (p are collected and all the others are combined into a total 
divergence. We may separate 

H = Ho + e(j) + Hp (53) 
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where Hq comprises all the terms not involving and then the polarisation piece Hp may 
be constructed from H — Hq — ecj). All gyroaveraging corrections (e.g., from 1 — Jq) are 
collected into Hp. We define the generalised vorticity VL and the polarisation vector P such 
that 



(54) 



sp 



The quantity on the right side is the gyrocenter charge density. Since the derivations are 
being done under strict quasineutrality (E'^/Sn) was neglected as discussed around Eq. [T3l). 
the quantity on the left side balances this. It is the polarisation density as developed by 
other methods in Ref. 28|. The generalised vorticity is defined in this manner as a quantity 



sensitive to small scales which (as seen below) in the MHD limit reduces to the simple ExB 
vorticity. 

Then the polarisation equation (Eq. [16]) is 

6fHp 

sp 



V p = -E 



(55) 



It is essential for the subsequent results to be able to write the polarisation equation in this 
form, with polarisation density given by a divergence. 

For any H with dependence on such that the separation H = Hq + e(j) + Hp yields 
dependence of Hp upon only through V0 and V]_4>, the functional derivative of fHp is 



5fHf 



/ 



dHf 



which can be written divergence 

5fHp 



dHr 



/ 



dHp\ 

dW<i)j 



dHp 

dW4> 



(56) 



(57) 



and then the species sum of the quantity in the square brackets is identified with P. Since 
within H only Hp depends on V0 or V5^0 we generally have 

dH ^ 



sp 



/ 



(58) 



and the need to be able to write Eq. fl55|l in that form is satisfied. The only restriction on 
the form of H is that all terms past first order in appear only through V0 or V5^0. Note, 
however, that there is no such restriction on A\\. 

For example, for the long-wavelength electrostatic H through second order. 



H 



2m 



mc 
252 



(59) 
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the polarisation equation is 
and the polarisation vector is 



sp ^ 



(60) 



(61) 



where we note that dW/B*^ commutes with spatial derivatives. In this case the species- 
summed velocity space integral is straightforward and we have 



n 



(62) 



where pM is the species sum of nm with n = J dWf the species density. One sees why the 
gyrocenter charge density plays the role of a (negative) vorticity, since V-P is proportional to 
— plus corrections due to the gradients of the densities and the magnetic field strength. 

The global conservation is trivial since the phase space integral of the gyrokinetic equation 
conserves particles for each species. The vorticity transport equation may be written as 



Here, the angle brackets denote the flux surface average, which is the same as the volume 
derivative of an integral over the volume enclosed by a particular flux surface 43|]. For 
simplicity we assume Hamada flux coordinates {V^C} where V is the volume enclosed by 



tions 



45|. 



the particular flux surface, the contravariant components of the magnetic field are func 
of V only, and the poloidal and toroidal angles (respectively) are unit-cycle periodic 
This leaves = 1. Then, = R ■ VV^ is the contravariant ^-component of R. 

We have used the property that the flux surface average of a phase space divergence 
annihilates the velocity coordinate derivatives (with respect to pz and /i) and commutes the 
integration / dW / B*^ past the spatial derivatives. Then, the flux surface average annihilates 
the angle derivatives in a flux coordinate representation, leaving the V^-component of the 
drift motion and the derivative d/dV . Note that the flux surface average of a kinetic quantity 
implies the species-summed velocity space integration, which is left understood. 

The introduction of the polarisation vector under the time derivative leaves this equation 
as a pure divergence. 



d dP^ 
dV\~dt 



+ feV 







(64) 
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where the superscript V denotes the contravariant V^-component. This describes charge 
conservation in the form V ■ J = 0, as a balance between the gyrocenter drift current (the 
species sum of /eR, including the parallel piece) and the divergence of P, so we may also 
identify dP/dt as the polarisation current. The quantity given by the flux surface average in 
Eq. which includes both pieces, can be taken to vanish everywhere given appropriate 
boundary conditions (e.g., it vanishes at the magnetic axis, where = 0, due to regularity 
of the vector component). Hence we will also have 

-,iV)[—^feVyo (65) 

where 7 is any flux surface quantity (also called flux function). Specifically, this expression 
vanishes for ■y = since is the quantity whose isosurfaces define flux surfaces, for a 
tokamak magnetic field. 



VI. A FURTHER COMMENT ON ORDERING 



In the conventional gyrokinetic ordering the small parameter is equivalently k\\/k± in the 
wavenumber anisotropy or the fluctuation amplitude e0/Te, which are used interchangeably 
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42l |. (The often stated ordering of ^ 1 is actually a posteriori as it follows from 
the ultimate requirement that the resulting dynamics is in the range of the diamagnetic 
frequency a;*, and for this to be small compared to the ion gyrofrequency requires p* <C 1.) 



Under these conditions we may simplify expressions by using the field aligned version 46|, 47 1 



of Hamada flux coordinates 



43 



-45|. This discussion follows the version used in Ref. [47 1 
which includes the definitions and construction algorithms for the coordinates. Starting with 
{l^^C} as above, we transform the toroidal angle only, defining = C — qd where q = q{V) 
is a flux function giving the ratio / in the contravariant components. Then, both 
B^ and B^ vanish, and the only nonvanishing component of B is B^ . This is defined as 
B^ = x' = dx/dV, where x = x(^) is another flux definition. The tokamak magnetic field 

B = /V(/? + VA^xV^ (66) 

may be written as 

B = VexVx = x'VexVV (67) 



which is called a Clebsch representation 46|]. The sign conventions are VRx'VZ-'Vip > and 



x' > 0. It follows that x' = —^ndip/dV and that for any vector (including the gradient) 
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the covariant components satisfy = = A^/2'k. With only nonvanishing, d/d9 
tracks the parallel derivative and therefore we may assume d/d6 <^ d/dV,d/dC, due to the 
wavenumber ordering. 

We consider the flux surface average of the toroidal momentum continuity equation (Eq. 

EOD 

where the focus is on the term in the drift motion involving A^/c, and Hf is the part of H 
which involves the fluctuations in the field variables. We may apply the flute mode ordering 
to the drifts term and assume B*^ = B + 0{p^,) dropping the small correction, so that 

^rr C F Y Cl OHf , . 

where we have used the exact equality F^^ / B"^ = which results from the index-raising 
operations on F^y = x' and the equality {B/x')"^ = 9^^ — 9^^ 9^^ giving the determinant 
of the perpendicular contravariant metric coefficients. Recall that = 1 and F^y = e^yeB^ . 
Substituting the drifts into Eq. (!68ll . the factors of (c/e) cancel and we find 

where the fiux quantities have been taken out of the fiux surface average. Finally we insert 
the 2iT normalisation in ^ so that 

We observe that the two fiux surface averages are the same term. Furthermore, for a local 
model we may evaluate A^ = A^y{V — Vq) near a zero at \^ = Vq because the amplitude of 
the fiux is arbitrary, and use 27iA^y = — x', so that 



_d_ 



{Vo-V){f 



OH/ 



where it is evident that the two terms are at the same order for any contribution by or 
A\\ to H. 

It is tempting to draw the cancellation obtained by differentiating the Vq — V factor. 
However, the remnant is of the same order. Though this manipulation does not yield a 
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useful final result it still serves to show that the potentially large term involving A^p|c gives 
a contribution at the same order as the terms we keep on the right side of Eq. (|50|) . Again, 
applying ordering to H only, the similarity in magnitude follows order by order, as the 
Hn contributions on each side are similar at each order. This dispenses with any notion 
that the large factor A^^jc should be accompanied by higher order drift terms due to if/ 
than are present in the term — (/if^,^) on the right side. We will find it necessary to keep 
the second order ExB energy term to obtain useful results, but since energetic consistency 
and the consistency in ordering we have just derived indicate, we do not need to keep terms 
beyond second order in and Ay on the left side of Eq. fISUl) . Again, the underlying energetic 
consistency yields straightforward conclusions consistent with the use of a single H (to any 
order) everywhere within any particular version of the gyrokinetic field theory model. 

The dependence of this result on energetic consistency cannot be overemphasised. In 
an ordering expansion, A^jc as a large term potentially introduces higher order terms in 
H. However, application of the field equations (through functional derivatives in ii) up 
to any particular order, which obey exact energetic consistency up to the same order, was 
used to arrive at the result that the two expressions involving {fH^^p) are at the same order 
(each contribution to if, order by order). Had the field equations been missing terms at 
the highest order kept in the drifts in the right side of Eq. (!50|) . this result would not have 
been obtained and we would be required to discuss spurious effects. This is why discussion 
of orderings in H past first order must be done in the context of energetic consistency. 

VII. THE TOROIDAL MOMENTUM TRANSPORT EQUATION 

We now do for toroidal momentum what we did for vorticity in Sec. |Vl Starting with Eq. 
(!50|) . we apply species summation and velocity space integration and then the flux surface 
average to obtain 

where as before V = H. ■ W, the velocity space integration annihilates velocity space 
derivatives and dW/B^^ commutes past the spatial divergence. We insert the definition of 
in Eq. (l34l) and collect the A^ terms to find 
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We pull out of the flux surface average to flnd 

Using Eqs. fl54f55p . the gyrocenter charge density is replaced by the vorticity under the time 
derivative, to flnd 

The divergence operator on P is done by parts to flnd 

+ |(feM + |7(/PAl') = -(/|^) (77) 

since V = W [d / dV) for any flux function. In pulling the divergence out of the flux surface 
average we have used (V ■ P) = d/dV (-P^)- We combine the flrst two divergence terms to 
flnd 

which vanishes by Eqs. dnHES])- Hence we have 

as a result for any H subject to the comment after Eq. flSHl) . through the dependence of the 
vorticity equation (Eq. [M]) upon the polarisation equation (Eq. [11]). Insertion of Eq. ( 17^ 
into Eq. (175|) and moving the right side term to the left side produces 

This result eliminates terms in which A^, the flux label, appears by itself (not under a 
gradient operator). We may identify the flrst and second terms as the time derivative of a 
total toroidal momentum density consisting of an ExB part and a parallel part (both have 
toroidal components). The next term gives the drift effects of the parallel part of the toroidal 
momentum density. The last term has already been seen to vanish under the total phase 
space integral, so as a flux surface average it should be another (set of) transport divergence 
term(s). It remains to show this. 
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A. Field Terms for Simple Hamiltonian Dependences 



We assume at first for illustration purposes an electrostatic case with general dependence 
of H upon and V0 irrespective of any ordering. We have the integrand in the last term 
of Eq. dZSD as 

(9^9 d(j) d(f dV(f) d(f 
noting that the partial derivative and the gradient commute. The gradient is done by parts 
to find 



dip \ d(f) dS/(j)^ 
The first expression in parentheses is the functional derivative and it vanishes under the 
eventual species-summed velocity space integration, due to Eq. flTB]) . The second term gives 
a divergence of a transport flux. Putting this back under the flux surface average, we have 

{/f>^l^{|--/^) (S3) 

which is the desired result. 

Following the same analysis for a general dependence of H upon and Ay, with the 
functional derivative with respect to A\\ not vanishing but replaced by the field term as per 
Eq. dn]) we find 

where the contribution due to dBj_/dip vanishes under the flux surface average. As we 
will see when explicitly showing the version of the result using the long-wavelength "MHD 
Hamiltonian" these two terms generally give the Reynolds and Maxwell stresses, respectively. 



B. Field Terms for Hamiltonian Dependences involving the Laplacian 

In the case that H depends on the field variables through not only and V0 but also 
V5^0 then there is more to do but generalisation is straightforward (for background see the 

48|). We expand 



text by Gelfand and Fomin 



dH ^ dHdl dH_ dl r_dl^y2 d± , . 
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again assuming the differential operator commutators vanish. The first two terms are done 
as before but the perpendicular Laplacian involves two integrations by parts. The Laplacian 
piece becomes 

dH 



9V]_0 dip dip 



Both of these pieces further expand according to 



and 



dip ' dV: 



-V 



'90 dH 
dip'^^^dWh 



dH 
dWc 



(86) 



(87) 



(88) 



The Laplacian piece is then 



dH 



r dH 2 ^ _ r 

■'dVl(l)^dip dip^^'dVl^ 



+ V- 



' dVl 



(89) 



Putting in the terms from the and V0 dependences, we have 

dH^dl(dH_ dIl JH_\ 
■' dip dip y d(j) ■' dV(f) dVl(j) i 



+ V 



/ 



dH 



9if 



(90) 



The first term in parentheses gives the functional derivative (which vanishes according to 
Eq. [T6l) and the rest become transport fluxes. The flux surface average with species-summed 
velocity space integration then yields 



90 / dH 



dH 



dip y dW<p ^^^•^9Vi0 ) ' ' ^ V-^*^ ' dW\<p^ 

which is the desired result. 

Following the same analysis for a general dependence of H upon and Ay, with the 
functional derivative with respect to A\\ not vanishing but replaced by the field term as per 
Eq. (dZD we find 

dH 



dH 



(91) 



90 / dH 

di^ r 9V0 



dip) ' dW] 



dv\ 

+ Vj 



dA 



dip 



' 1 dH dH 



'Mil ^ dH 



dip ' dV\A\\ 



(92) 
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where the contribution due to dBj_/d!f vanishes under the flux surface average. 

Though more comphcated, this result shows that for general dependence of H upon </> 
and A\\ and their gradients and Laplacians, all terms in the toroidal momentum transport 
equation may be recast as divergences of flux surface-averaged transport fluxes, and no terms 
with not under gradients appear. 



C. Illustration using an MHD Hamiltonian 

We may dispense with finite gyroradius effects while studying equilibrium flow dynamics 
in the long wavelength regime (typically for edge equilibrium flows the local is between 
1/100 and 1/30) and also expecting > (l/nee)V^Pi unless the flows are weak I27I. In 
conventional tokamaks the low-frequency and low-beta assumptions referred to above (Eq. 



[T^ are wel 
(cf. Refs. [7 



satisfied and the relevant form of MHD is reduced MHD 

n. 
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35|. In general 



lll| ) there are contributions to polarisation by but in this limit the terms 
due to A\\ in are small compared to B\/8tt, so for these purposes it is sufficient to 
keep A\\ only in the parallel kinetic energy term mU'^/2 and in B^/Stt. These considerations 
are interesting in their own right and will be treated in a different work. It is presently 
more important, however, to concentrate on the structure of the theory rather than the 
details of any particular version, so here we use the one necessary to obtain familiar reduced 
MHD forms. One further consideration deserves emphasis: no equation beyond the original 
statement of L and H requires justification; only L and H themselves. Once L and H are 



chosen, all assumption stops and 



The Hamiltonian used in Ref. 



urther results are a matter of derivation. 



• 3 



IS 

f/' „ 1 



H = m h /i-B + e0 mv% (93) 

where the square of the ExB velocity and the parallel velocity are 

4 = ;^|Vx0r mU = p.--A, — = U (94) 

For the derivatives of H we have 

dH dH mc^ dH e dH 

Vx0 7r^ = —U 7^ = (95) 



dcj) dW(t) 52 dA\\ c dWA 
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and there is no dependence on the Laplacian of either of the field potentials. The functional 
derivatives are 



6fH 



fe 



The polarisation vector is 



1 „ ^ fmc^„ 
— V ■ Bu - V 



Mil 



- j dW-JU (96) 



(97) 



Inserting these forms into Eq. (IHOjl using Eq. ( 18^ . we find upon collecting the W terms 
d / c 



dt 







(98) 



which is the desired result. We have the conserved local toroidal momentum density, the 
quantity under the time derivative. The pieces are the covariant toroidal angle components 
of the ExB and parallel momenta, respectively. The quantity under the divergence is the 
radial flux of the toroidal momentum. The pieces are the magnetic flutter and ExB/parallel 
Reynolds stress, the pure ExB Reynolds stress, and the Maxwell stress. This is the result 
as obtained without the use of any ordering except that involved in the prescription of L. 

We also consider this as a mean field theory, in which the flux surface average is also 
understood to contain a time average over a mesoscale range longer than eddy correlation 
times but shorter than transport diffusion times. The fluctuations are assumed to be small 
amplitude (order in relative amplitude, with velocities normalised to the sound speed). 
Following the considerations in Sec. IVTI we apply flute mode ordering to the fluctuations, 
with the observed relationship between ip and ^ in covariant components. Turbulent fluxes 
appear in the flux surface averages where the ^-component of equilibrium flows is negligi- 
ble (the magnetic drifts constitute neoclassical transport, which we neglect here). In the 
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for 



transported quantity the contribution due to the fluctuations is neglected. See Ref. 
application of this in a drift-fluid model. We will also take the single-fluid MHD approxima- 
tions, where each species is assumed to have the same parallel velocity and in the end there 
is a single pressure. In this section, the tilde symbol denotes fluctuations in the indicated 
quantities. 

To see the ExB toroidal momentum we note that 



fi2 



(99) 
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using Eq. dnnD and R^\V^\ = 1 with R the toroidal major radius, and noting that the 
toroidal magnetic field does not contribute. Thus we may combine 

Pm-|jV0 ■ VA^ + fp,b^^ = (pmu^) (100) 

into a toroidal momentum density given by the mass density times the toroidal fiow u,^ 
under the mean field and single-fiuid MHD approximations. 
For the Reynolds and Maxwell stresses, we observe that 

^X7V V.I- ^ 9'''' d4> g""^ 

Since F^^ / B'^ = the ExB velocity components are 

~v c 90 ~£ c dcf) 

X X dV 

The covariant metric coefficients are 

gVV _ _gV^ 

^« = gVVgi^_gVigiV ^ = ^yVg^^ _ gV^g^V ^^^^^ 

Using the index lowering operation 

= g^i^^ + g^vv^ (104) 

we have 

-;J||V\/-V0 = t5^ty^ (105) 
and putting the 27r normalisation back in, we have 

-;^^vv-v0 = ?y% (106) 

which is the radial/toroidal Reynolds stress. With magnetic ffuctuations given by 

= -1^ = (107) 

we similarly have 

^^VV ■ Vi|| = -^6% = -5^5, (108) 
which is the radial/toroidal Maxwell stress. 
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For the magnetic flutter nonlinearity and the ExB/parallel Reynolds stress we note that 
the part of the drifts due to fluctuations is given by 

~ cF ~ cF ~ F 

R = Vi/— = V0--^-f/VA||.- (109) 

where to lowest order in we approximate — i- i? and neglect the second order drift term 
due to —mv'^/2 in H. The \^-component of this is 

The flux surface average of the fluctuation drifts term is then 

{fP.Kv) ^ {fpX^^ - fUpXf.^) (Ill) 

These are two different effects in a gyrofluid sense because they involve different moments 
of /. In both terms / — )■ / represents fluctuations because a flrst-order fluctuation term 
vanishes under the time average contained in (). Considering this, we have the moments 

(nmu\\v^b^ + p\\b^b^'^ (112) 

where and are as given in Eqs. fll02fl07p . These two terms give the ExB/parallel 
Reynolds stress and the magnetic flutter transport pieces. The distinctions between U and 
Pz/m and between and Py = + nmuj are neglected for small amplitude fluctuations. 

At the MHD level the sum over species leads back to the total mass density pm and the 
total pressure p which replace nm and p\\ {i.e., neglecting anisotropy), respectively. This 
leads to 

g d I ^ ^ 1~~_~\ 

^ (Pmu^) + -g^ (pmv^v-o + pMU\\v^b^ - —B^B^ +pb^b^) = (113) 

under MHD and mean fleld approximations having used the "MHD Hamiltonian" given in 
Eq. fl93|) . We have placed the two Reynolds stress terms next to each other and left them 
separate, for clarity. All of these terms are easily identiflable with well known processes 
within the MHD fluid model. We note that the mean fleld ordering has been used for 
evaluation purposes only, with the actual transport equation for this model given by Eq. 

This exercise has served to prove that the general gyrokinetic toroidal momentum con- 
servation laws can be brought back to MHD via straightforward application of the MHD 
approximations. 
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VIII. THE ENERGY TRANSPORT EQUATION 



The energy transport equation is much better known, even for gyrokinetic theory. How- 
ever there is an important manipulation concerning the part of H due to 0. Decomposing 
H = Hq + e(j) + Hp as in Sec. |Vl we start with 



fH = fHo + fe(l) + fHp 



Summing over species and inserting Eqs. fl54|55p we find 



;ii4) 



sp sp 



;ii5) 



and under the flux surface average 



{fH) = {f{Ho + Hp))-{P-V<P) + 



_d_ 

dV 



V 



(116) 



We also have 

{fHV) = {[f{Ho + Hp)-P- V0] V) + ((V ■ 0P) V) (117) 

where as before = R ■ W. Reduction of the time derivative term in Eq. f lSTj) is done the 
same way as in Sec. IVllI for the toroidal angle derivative term in Eq. ( l50l) . noting that the 
time derivative of Bj_ survives, as in the analysis leading to Eq. (l33l) . We find 

d 



d_ 

di 



'fHo + fHp~{P-V<f)) + ^\ + 



dV 



fHoV 



^(ifHp-{P-V^) + V-^P)V + ^^ 



V 



dv\ 
dV\ 



(9V 



dH 



dt ' 9Vi 
dH \ 



dt V47r 



dV^ All 



+ Vj 
= 



'Mil ^ dH 



dt 'SV^Aii 



;il8) 



In practical cases only the terms on the top line are significant. Under the time derivative we 
have the thermal and kinetic energy, ExB energy, and magnetic energy. The transport term 
gives the ExB advection (the V0 term in H in the drifts) and magnetic flutter (the f/VAy 
term in H in the drifts). All the others are polarisation and induction corrections. Inserting 
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the MHD Lagrangian from Sec. IVII C\ and taking the single-fluid MHD approximations, we 
flnd 

{fHo) = 1{P) + 1 (PMUI) {fHp - P . V0) = ^ (pMvl) (119) 

together with B'j_/87c for the energy pieces, and 

(fHoV) = ^ {pv"") + (120) 

where ge|| is the conductive electron parallel heat flux, for the dominant transport pieces. 
The correspondence with MHD forms is evident. This exercise has served to prove that the 
general gyrokinetic energy conservation laws can be brought back to MHD via straightfor- 
ward application of the MHD approximations. 



IX. SUMMARY AND DISCUSSION 



The main results of this work are the antisymmetric 4-bracket form of the gyrokinetic 
equation, Eq. flM|) . the global conservation laws for energy and toroidal momentum, Eqs. 
fl33|48l) . the local phase space advection equations for energy and toroidal momentum, Eqs. 
fl28|4ip . the phase space continuity equations for energy and toroidal momentum, Eqs. 
(15011511) . and the transport equations for vorticity, toroidal momentum, and energy, Eqs. 
(I63f80llll8p with Eqs. ( I8^f92!) as auxiliaries. It is important to note that no ordering as- 
sumptions were required to obtain these. The only required conditions are the form of the 
Lagrangian and Hamiltonian implied by Eqs. (fT|2l[T3|) . the ability to write the vorticity in 
the form of a polarisation vector divergence as in Eqs. (I54|55p . which enable the cancellation 
given in Eq. (1 79 p . This in turn leads to the single condition on the form of that the 
dependence upon must be as in Eq. ( l53l) . which leads to Eqs. ( I5^55ll58p . These results 



confirm what Refs. 



12 



13| already implied and Ref. ll|] already reviewed. The part of the 
results that are novel comprises the 4-bracket form of the gyrokinetic equation, its use in 
proving the conservation laws, and that the local forms of the conservation laws, i.e., the 
transport equations, have a solid fundamental basis for any reasonable choice of Lagrangian. 
The correspondence to nonlinear reduced MHD was shown, also without ordering, in Eq. 

(EHD. 

The usefulness and significance of the ability to arrange the Lie transforms with which Lp 
and H are built such that all time and toroidal angle dependence is kept in the dynamical 
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field variables and the latter are strictly contained in H {i.e., to obtain the form of Eq. [T]). 
Previous work has already shown how to arrange this, for both the conventional small scale 
and also for the newer large scale orderings js, 9|. 

The role of the cancellation in Eq. f l7^ highlights the role of the vorticity equation in 
the overall consideration of momentum. In a low-frequency fluid drift model, the natural 
decomposition of rotation is not poloidal/toroidal but perp/parallel, with the ExB energy 
equation controlling the evolution of the perpendicular flow. In a gyrokinetic or gyrofiuid 
model, this role is taken over by the ion gyrocenter density variable, the energy content is 
controlled by JHe (where He is the part of H which depends on 0), and the relevant con- 
served quantity is, one and the same, gyrocenter (apparent) charge density and generalised 
vorticity. The appearance of the polarisation density as the divergence of a polarisation 
vector is a fundamental property underlying the conservation laws. Hence, the local con- 
servation of toroidal momentum depends on the simultaneous conservation of vorticity. At 
the level of the equations, obtaining the one is dependent on the use of the other to pro- 
vide a cancellation by which the amplitude of the poloidal magnetic flux (as opposed to its 
gradient) is removed from influence on the transport of toroidal momentum. 

The mathematical functioning of the global conservation equations for energy and toroidal 
momentum depends critically on the ability to employ time and toroidal angle translation 
symmetry, respectively, in the use of functional derivatives to evaluate an apparent residual 
which in fact vanishes under the integral (the right hand sides of Eqs. [281I4T]) . This in turn 
depends on the exact working of the Euler-Lagrange equations for the field variables, and 
this in turn is defined in terms of functional derivatives. 

At the level of the local equations, these apparent residuals are shown to be recastable 
in terms of divergences using techniques similar to the use of functional derivatives. Hence, 
for any H, they are recast as transport fluxes. All of the transport fluxes which remain 
significant under mean field and small fluctuation ordering are found to have clear MHD 
analogues, another step in confirming the general solidity of the theory. Using a simplified 
Hamiltonian defined in Eq. ( l93l) which keeps finite gyroradius residual corrections only in the 
retention of the ExB energy as a second order drift, it was shown to be straightforward to 
recover the MHD limit of the transport equations (Eq. [98l) . with the familiar Reynolds and 
Maxwell stresses and the magnetic flutter nonlinearity elucidated by a mean field analysis 
(Eq. I113p . This "MHD Hamiltonian" provides the link between gyrokinetics and nonlinear 
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reduced MHD at the total-f level. The correspondence at the delta-f level was recently 



proved elsewhere for gyrofiuid theory 



49|. 



A subsidiary analysis was used to show that for small scale fluctuation ordering, the 



magnetic flux term (A^/c), which is formally of order e^^ as discussed in Ref. jSj does not 
introduce terms in H at higher order than already are necessary to evaluate the momentum 
transport equation. This result, in Eq. (172|) . is sufficient to allay recent concerns about the 



integrity of the treatment o: 
have been voiced by others 



momentum conservation and transport by gyrokinetics which 



3l|. 



Second order drifts are necessary for energetic consistency if polarisation is nonlinear, 
but it is possible to have a model with linearised polarisation and only first order drifts 
in the gyrokinetic equation itself and still satisfy energetic consistency. Such models also 
have well behaved local transport equations for toroidal momentum, as shown in Appendix 
IB 1[ The second order Hamiltonian terms — quadratic in the field dependent variables — 
must be present in some form, either as second order drifts in the gyrokinetic equation or 
as background field terms, since they control the dynamical energy accounted for by those 
variables. 

Energetic consistency underlies all the derivations, which do not work otherwise. Any 
model which does not have this at its heart and which attempts to go beyond first order 
field dependence in drift terms is prescribed to fail on consistency grounds. It is important 
to note that the abovementioned concerns were not done with an energetically consistent 
analysis and did not address any of the previous results on energetic consistency nor even the 
existence of gyrokinetic field theory. However, the seminal references on gyrokinetic field 



theory already clearly demonstrated the need to keep fully nonlinear po 
level of ordering past first order dependence of H upon the field variables 



arisation at any 



13|, and this 



was and is followed rigorously in field theory treatments of reduced (drift) fiuid equations 



5G 



Gyro 
models 



511] and total-f gyrofiuid equations |36l. l37|. 

dnetic field theory was not necessary to build the original gyrokinetic computational 
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301 1 . However, an equivalent field theory can be built from the given gyrokinetic 



Poisson equation (multiply by the time derivative of (f) and find the resulting form with / 
times the time derivative of H, then use that H to re-derive the equations). The above-cited 



computations and more recent ones derived from them 
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456( 1 are consistent with this. This 



is why conventional delta-f global computation is on solid fundamentals. 
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Gyrokinetic field theory is however the only form of modern gyrokinetic theory as it 
is required to provide well founded generalisations in practically any context. It is also 
necessary to any effort which aims to generalise gyrokinetic theory in practically any way. 
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Appendix A: Use of Functional Derivatives 

The Euler-Lagrange equations for the fields are found by identifying the functional deriva- 
tives of L with respect to each field potential and setting them to zero. For an introduction 
to the use of functional derivatives in this type of field theory see Section II B of the review 
by Morrison js^]. For background on functional derivatives within the general context of 



the calculus of variations see the text by Gelfand and Fomin [48|. For the more familiar 
version by which stress tensor theorems are proved within classical field theory see the text 
by Landau and Lifshitz 3^. 

The functional derivatives are defined with respect to the space covered by dV, so that, 
e.g., variation of results in 

5L{5(t)) = I dV5(t)— (Al) 



in which it is understood that the definition of the functional derivative includes velocity 
space integration. The arbitrariness of 5(f) and the extremal requirement for L give the field 
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equation as 

Since the dependence of L on is through H only then this is equivalent to 

E^ = (A3) 

Variation of A\\ results in 

5L(M||) = j dVSAii (A4) 
and the arbitrariness of 5A\\ and the extremal requirement for L give the field equation as 

Since the dependence of L on A\\ includes the field term —B^/Stt this is equivalent to 

op II 

These are the field equations for arbitrary H given the field term B^/Sn. The induction 
equation depends on the form of this term; for example the version of Bj^ which gives a 
result closer to the form of Ampere's law used in the Grad-Shafranov equation for MHD 
equilibrium is 

^i = ;^|v±(A||i?)f (A7) 

where R is the toroidal major radius. In this case the field equation is given by 

op II 

with the new form on the right side. Due to the axisymmetry of R, either of these forms 
may be used in the above derivations. 



1. Integrals over / and H 

The conservation laws depend on phase space integrals over expressions of / and H in 
the form 

Y^jdKfdaH (A9) 

sp 
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where da is a differential operator with the hnearity property, (e.g., a time derivative). In 
turn, the dependence of H on the fields involves the field amplitudes and their derivatives 
which we can schematically express as 

H = H{(j), db(j)) (AlO) 

where dh is another differential operator with either the linearity or the the Hermitian prop- 
erty (e.g., V or V5_ or Jq). The relevant differential of H is 

OH dH 

daH = -Tr-da(t) + ^TT (All) 
ocp odh(p 

We can turn the integral into a functional derivative form if and only if the two operators 
commute: 

JdAfdaH = JdA ^ {dadt - d,da)<P = (Al2) 

If all the (if, in H commute with d/dt then we have time symmetry which is a prerequisite 
for energy conservation. If all the di, in H commute with 8/ dip then we have toroidal an- 
gle (axi-) symmetry which is a prerequisite for toroidal momentum conservation. In general 
the commutator does not vanish if da is V, since db involves B which is spatially depen- 
dent. In stellarator geometry dB/dip 7^ and therefore toroidal momentum conservation 
is not conserved (it is exchanged with B). In tokamak geometry poloidal momentum is 
exchanged with B but toroidal momentum is conserved in the ideal case (no destruction of 
axisymmetry) . 



Appendix B: Equivalent Hamiltonians to Conventional Models 



As an example of how one can back-construct a gyrokinetic field theory from a con- 
ventional model, we consider the gyrokinetic Poisson equation as currently used in most 
numerical simulations 
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54 



56| . The gyrocenter charge density piece is kept with gyroav- 



eraging and the polarisation is given by a background ion density 



^ / dW[eJof] + V 

sp 



B^ 







(Bl) 



For this to come from a Lagrangian it has to be the functional derivative with respect to 
We may rebuild this by multiplying by 5(j) and integrating over the volume, so that 



Y^f dkdct) [e Jo/] + [dV6(t>\/- ^^^^ Vx^ 



sp 
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(B2) 
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Here, all we need know about Jo is that it is Hermitian. Doing the relevant integrations by 
parts we have 

Y^fdA [feJom] - fdV ^^^Vx0 ■ Vx(50) = (B3) 



sp 
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We identify the relevant parts of the Lagrangian as 
L = 



sp 



252 



(B4) 



Variation of this with respect to indeed recovers Eq. fIBip . Inspection of the version of the 
gyrokinetic equation used in these references shows that indeed He = eJo0 is used as the 
perturbed Hamiltonian in the drifts. Therefore, energetic consistency is assured. 

Most versions of this model are considered with electrons in adiabatic force balance 
parallel to the magnetic field. 



dWfe = no 



(B5) 



with gyroaveraging neglected and with the same no as in the background ExB energy term. 
The factor cq is the fundamental unit of charge, so written to distinguish from use of e 
as a species charge. The temperature Tg is a flux function and the subtraction of the flux 
surface average reflects the vanishing of the parallel gradient for this component. This state 
of force balance is expected to evolve adiabatically in local thermodynamic equilibrium, so 
the contribution to L is a profile anchor piece plus a fluctuation energy piece. 



-f'electrons = noeo4> + ^o^ {(p - 



(B6) 



which becomes a field term. With this substitution made in L, the polarisation equation is 



no— {(f) - (0)) - V ;;t^V_ 



-noeo + J2 / dWeJof 



(B7) 



'Te '"^ ''''' ' 

where the sum is over the ions only (for singly charged ions the same no is used everywhere; 
for more species the appropriate adjustments are made to keep the profiles charge neutral). 
The term noCo subtracts the profile piece from the ions. Energy conservation for this model 
was proved for the ORB code in Ref. |5J]. The same model is also used by the GYSELA 
code 



The form given by Lee in the original gyrokinetic/Poisson system 28|, l29|, ISSj is 



2TQ^ 



50(R) 



OR 



(B8) 
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where \l/ is the potential with which particles are pushed in a gyrokinetic particle in cell 
model, (f) is given strictly as a function of the gyrocenter position R, the symbol used for 
gyroaveraging is replaced by our use of Jq, and q, T, and (ft/fijp = mTc^ jq^B"^ are constant 
parameters, with q the same as our e. This is Eq. (2) of Ref. [29]. It is also an electrostatic 
model with A\\ = and Pz/ttl interchangeable with a parallel velocity. 
In our notation we write this as He = e\l/ with 

2 

HE = eJo<P-^\V<pf (B9) 

In other words, this is the same choice as our MHD Hamiltonian (Eq. [93]) . except for the 
use of the gyroaveraged as the first term and the terms dependent upon A\\ in the MHD 
case. Lee's gyrokinetic Poisson equation is written in Eq. (3) of Ref. [29| as 

2 

- T^(l - To)(p + ■ = -47reo(n, - n,) (BIO) 

Ad Xjj riQ 

The operator Tq results from Jq integrated against a Maxwellian and is given formally by 
multiplication in wavenumber space of Fourier coefficients by the Tq^B) = e~^/o(6), where 
Jo is the zeroth modified Bessel function, and the argument is 6 = k'j_p'^ with p the species 
thermal gyroradius given by = mTc^/e^B^. The factor rTj with the overbar denotes 
velocity space integration of Jofi, the rii without the overbar is velocity space integration 
of fi, and t/X^) = AirnQel/Ti and p^Z-^D — 4:'7inQMiC^ / Bq and no are constant parameters. 
Using the definitions in his Eqs. (4-7), this is found to be equivalent to 



(Bll) 



where is a Maxwellian with species parameters n and T, the correction factor including 
Jq and k^p"^ restores full finite gyroradius (FLR) effects to the field term, the sum over 
species using mg — )■ includes the electrons only in the last term e Jq/, and we have restored a 
toroidal model with the factors B and B*^ arranged to preserve Hermicity of all the operators. 
The first term is true charge separation, the second with F*^ is a field term, and the rest is 
the part due to the dependent variable. If we prescribe an electrostatic Lagrangian, 

" ' mc'^ ,„ ,,n] /■ „ . 1 



L = E/^A/^. + E/^AF^^ 1^ [02 - (Jo0)2] - ^ |V^0r| + / dV^rn' 

(B12) 



sp sp 
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with the particle Lagrangian given by 

Ln 



e \ 2m 



(B13) 



then the Euler-Lagrange equations recover gyrokinetic equations with e\l/ as the 0-dependent 
part of H and the polarisation equation given by Lee as his gyrokinetic Poisson equation; 
in other words, Eqs. (1-3) of Ref. |2^]. These were given for slab geometry but with the 
dependences of B placed as shown here together with the use of Bq where indicated, this 
recovers a toroidal model with some extra field terms in L compared to ours. Since the field 
correction term in the polarisation equation is the functional derivative of a positive definite 
quantity, it still conserves energy. Since we are able to recast the field terms as a field 
Lagrangian, and to show otherwise that the particle pushing potential and the gyrocenter 
charge terms arise from the same term He in the Hamiltonian, energetic consistency is 
entirely satisfied. 



Various simplified versions are given. For example, in the GTC code 59|, |60[ the particle 
pushing potential is consistent with He = eJ^cj) and the gyrokinetic Poisson equation is 
given as 



;i - ro)0 = 47reo(Tii - n^) 
then we have a linearised polarisation model consistent with 



\2 ' 



sp 



L = Y.ldAfL, + Y:Jd\F''^ 



(B14) 



(B15) 



where Lp is as in Ref. (IBlSp with the simplified He^ and adiabatic electrons are prescribed 
through Eq. (IB5|) above. Since using this L the above polarisation equation is recovered along 
with use of ^ = J^cf) in the gyrocenter drifts, the model remains energetically consistent. 
Approximation of the F^'^{Jq — 1) to long wavelength then recovers the ORB/GYSELA 
model mentioned above. 

Lee mentions these models also in Ref. 



28 



giving their source in another model used to 



obtain his Eqs. (21,22). The particle pushing potential is given as 



^^/ = Jo(j)+ — 



M 



(B16) 



in our notation. The factor 1/T is obtained by approximating df /dfi as {—B/T)f since the 
second order term leads to the polarisability, which should be dominated by the largest scale 
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part of / which should not depart significantly from a Maxwellian. This \& was then used 
to recover the model given in Eqs. fIBSIIBlOl) by using the long wavelength approximation of 
the second order term. We write the interaction part, JHe = /e\l/, of the Lagrangian as 



Y^jdKfHE (B17) 



L = 

sp 

then undo this effective integration by parts under the phase space integration to obtain a 
similar model with 

He = eM - [M4>') - iMf] (B18) 

With this form the field term FLR correction is unnecessary and with quasineutrality we 
recover use of an L in which the only dependence upon is in H. The corresponding 
polarisation equation is 

Y^JdW [eJof + {JoMJo - [JoM]) 0] = (B19) 



sp 



where the polarisability Ai is given by 

df 

M ^ --± (B20) 

It is very difficult to be able to take this derivative in a particle in cell model, which is why 
this second order term is usually approximated in some fashion. 



This is very close to the version given by Hahm 6| , who obtains an additional term from 
the second order part of the Lie transform so that 

= e (0) R - 1^1^ (0') R - ^ (V$ ■ bx V0) R (B21) 
where in this case the subscripted angle brackets denote Jo at the particle level and 

= 0-(0)r ^ = j'^d'd (B22) 

The first of the two second order terms is equivalent to Lee's. The second involves an indef- 
inite gyrophase integral and is difficult to compute. The equivalent polarisation equation is 
then given keeping all the relevant terms through use of the Lie back-transform. As noted 



above, this is equivalent to the use of field theory However, the added term in He 

yields a contribution to the polarisation equation which is one order down in from the 
others, and it is never kept in computations. Nevertheless, with the approximations always 
effectively made in L and nowhere thereafter, exact energetic consistency is preserved in all 
of the versions. This effectively brings us back to Lee's He- 
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1. Momentum conservation with linearised polarisation 



Here we demonstrate the existence of a familiar toroidal momentum conservation equa- 
tion within the simplest possible form of a global gyrokinetic model. The dynamics is 
electrostatic, quasineutral, gyroaveraging is neglected, and polarisation is provided by a 
background field term. The Lagrangian is 

L = Y.jdKfL, + JdV^lV^^pf (B23) 

sp 

where pm is given by a static profile (or constant), and Lp is of the form given in Eq. ([1]), 
with 

2 

H=^ + fiB + e(p (B24) 
2m 

The corresponding polarisation equation is 

V ■ ^V^0 + Y^JdWef = (B25) 

This is the same model as in the ORB code mentioned above, except we simplify by neglect- 
ing (1 — Jo) hence FLR effects. 

The particle drifts are given by Eq. (j4]), with spatial part 



Bm = \/^-^ + fiVB-~-^(v-~]+ (B26) 



with the pieces identified as the ExB velocity, the grad-B and curvature drifts, and the 
parallel velocity. The radial component of this is then 

V = ^VR= {veV + {v^bV + {vcV (B27) 

which we use below (the parallel piece does not contribute). 

The steps to the momentum equation are the derivations of Eqs. fl41ll50ll63p which remain 
unchanged, and the polarisation vector 

- P''''\^<P (B28) 
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with the only difference to the MHD model at this point being the static pm- Specifically, 
we still have 
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with the background pu- The cancellation in Eq. (179|) remains intact, and we have Eq. (IHOj) . 
We evaluate 

to obtain 

^ {Pm{ve)^ + M^) + ^ ( ■ (^/p.&^R - Pm;J|^V0^ ^ = (B31) 

as the local toroidal momentum conservation equation, i.e., the toroidal momentum trans- 
port equation. The pieces are easily identified as the ExB and parallel components of the 
toroidal momentum, the ExB/parallel Reynolds stress, and the toroidal component of the 
pure ExB Reynolds stress, all familiar effects. In the drifts, {ve)^ gives the ExB/parallel 
Reynolds stress, while (fys)^ and (fc)^ give the grad-B and curvature drifts which comprise 
the neoclassical transport (to evaluate these, of course, one needs treatment of the collisions; 
cf. the review by Hinton and Hazeltine j43|). Except for the corrections arising through Jq, 
this is the toroidal momentum equation satisfied by the model used in the conventional 
numerical models mentioned above. 

The recent NEMORB code [61;|] is the electromagnetic version of ORB. It reestablishes 
A\\ at first order in H, placing the second order term alongside the polarisation field term. 
The version neglecting FLR effects is given by 

i - E / / + / I V.,r ^g;rA-i I V.A„ r) (B32) 

with reduced mass M = rrieMi/ {nif. + Mj) and with Lp as in Eq. ([1]), and H given by 

H = ^ + liB + e((t)-^A^ (B33) 
Im \ mc J 



The polarisation equation is unchanged, given by Eq. (1B25|) . and the induction equation is 



These are the simplifications of Eqs. (14,25) of Ref . respectively, descending accordingly 
from the simplification of H (see the discussion at the beginning of Sec. IVII Cp . 
The corresponding toroidal momentum transport equation is 

d d I f 86 1 dAw \ \ 

^ {Pm{ve), + fvA) + ^ ( • [fpzb,K - PM^^ V0 + ^^VA|| j ^ = (B35) 
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which merely adds the Maxwell stress as the last term. Other terms from the field equations 
are put into the form dS/dip with scalar 5* and are annihilated by the fiux surface average. 
The neglect of FLR effects in this demonstration is merely for clarity; we may restore Jq and 
expand it in a series of Laplacians to recover several correction effects following the procedure 
for Laplacian field variable dependence given in Sec. IVII B[ Hence, we have demonstrated 
that momentum conservation and transport in ORB and NEMORB, as well as other related 
codes, is on a solid foundation. 



Appendix C: Torque due to a Charge Source 



In exceptional cases the particle sources are not in charge balance. The exemplary case 



is ion orbit loss 
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63| . in which the ions on large trapped (banana) orbits impact material 



surfaces while the electrons remain confined. Strictly speaking this is a transport effect and is 
accounted for by the drifts term and a loss flux through the boundary. In practice, however, 
it is modelled by a localised loss term in the ion continuity equation. The corresponding 
vorticity transport equation is 

dfi 



Sc 



dt 



(CI) 



Since Sq is a scalar quantity we may specify 



(C2) 



with s = ds/dV = on the magnetic axis {V = 0), without loss of generality. With 
V ■ P = —Q we have 

f) I r)P^ . \ 





d I dP^ ■ ^ 



SO that the flux surface average quantity vanishes as before. 

In the toroidal momentum continuity equation the terms which involve are 

1. 



The terms on the last line vanish. With the manipulations which flnd 



-Vs ■ VA 



(C3) 



(C4) 



(C5) 
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we find 

It <^'>+l^t (^'^)+t - m (-^^'^x'' ■ V*')-(i«^BxV. . V,) (C6) 

In the case of tlie MHD Hamiltonian P = —Pm{c'^/B^)'V±4' the first term gives tlie 
ExB covariant toroidal momentum (equivalent to toroidal angular momentum). Hence the 
second term can be represented as a toroidal torque 

T=-^BxVs R^T-Vip^T^ (C7) 

and the vorticity source is given by 

5o = V-^TxB (C8) 



47 



